Decoherence of charge qubit coupled to interacting background charges 



Igor V. Yurkevich,^' ^ Jim Baldwin,^ Igor V. Lerner,^^'^ and Boris L. Altshuler'*' ^ 

^School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, United Kingdom 
® Institut fur Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 
"'institute for Nuclear Theory, University of Washington, Seattle, WA 98195, USA 
^Physics Department, Columbia University, 538 West 120th Street, New York, N.Y. 10027, USA 

The major contribution to decoherence of a double quantum dot or a Josephson junction charge 
qubit comes from the electrostatic coupling to fluctuating background charges hybridized with the 
conduction electrons in the reservoir. However, estimations according to previously developed the- 
ories show that finding a sufficient number of effective fluctuators in a realistic experimental layout 
is quite improbable. We show that this paradox is resolved by allowing for a short-range Coulomb 
interaction of the fluctuators with the electrons in the reservoir. This dramatically enhances both 
the number of effective fluctuators and their contribution to decoherence, resulting in the most 
dangerous decoherence mechanism for charge qubits. 
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An implementation of quantum logical gates using 
solid-state nano-devices looks rather desirable due to 
their potential scalability. One of possible routes in this 
direction is developing of charge qubits. Schematically, 
a charge qubit is a contact of two normal or supercon- 
ducting islands with charge carriers tunneling between 
them. The gate-controlled Coulomb blockade separates 
two charge states from all the others thus making a qubit 
P) 0) S| ■ Experimentally such charge qubits have been 
built as double quantum dots (DQD) 0, H, 0, 0] or 
Josephson junction structures, made of a small super- 
conducting island connected via a Josephson junction to 
superconducting reservoirs 
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Decoherence due to coupling to the environment re- 
mains a major impediment for developing charge qubits. 
This determines the demand for a better theoretical un- 
derstanding of the microscopic mechanisms that lead to 
decoherence. It is widely believed that the main contri- 
bution to decoherence in charge qubits comes from its 
coupling to fluctuating background charges (FBCs) that 
create dynamical electric field, which affects qubit charge 
states (see ref. for reviews of experimental evidence) . 
Since charge impurities are spatially quenched in the ex- 
perimental temperature range T ~ 20 ^ 50mK, the most 
probable cause for FBCs is from random recharging of 
electron traps embedded in an insulating layer close to 
the electronic bath. The role of the bath can be played by 
normal as well as superconducting metallic leads or gates 
when the electron states on the traps are hybridized with 
those in the bath. The appropriate model for FBCs has 
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been developed [lil [igI [itI 

analogy with the spin-fiuctuator model of the spectral 
diffusion in glasses [24|. In its frame, relevant experimen- 



tal findings on the decoherence and dephasing in charge 
qubits have been successfully explained [25l |. 

The question which has not yet been addressed is 
whether there is enough of FBCs to cause experimen- 
tally observed decoherence. For a typical DQD qubit, the 



area of the electrodes is of order of firri^ and the volume 
available for the FBCs which can contribute to decoher- 
ence can be estimated as 10^^^ 10~"'^^cm'^ so that their 
number is A'goom ~ (10 100)c where c is their number 
concentration measured in ppm (parts per million). 

As the appropriate electronic levels are randomly 
spread in the energy interval of the order of the band- 
width, D ^ 1^ lOeV, the effective number of traps iVoff is 
reduced in comparison to iVgoom by a factor of S / D where 
5 is the interval of energies in which FBCs are effective. 
Classical considerations l^, 3, 13] suggest that 5 is of 
the order of the relevant temperature, T ^ 10 50/iK, 
while the electronic trap levels separated from the Fermi 
sea by energies £q ^ T are frozen out, i.e. their con- 
tribution is exponentially suppressed. Had this been the 
case, would be five orders in magnitude below A^goom- 
Then the probability to find even a single effective FBC 
would be negligible for any realistic value of c. 

23| show that 
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The full quantum considerations 
electron traps energy levels up to £q < 70, contribute 
to decoherence. Here 70 is the hybridization energy, i.e. 
the broadening of the appropriate level. Its typical value 
can be much higher than T but still is likely to be much 
smaller than D. The effective number of traps is thus 



(1) 



Even with this enhancement N^s still seems to be of or- 
der 1 or even smaller for cleaner samples, with c ~ Ippm. 
At the same time experiments [8-13] show that even the 
cleanest samples suffer from considerable decoherence. 
Moreover, assuming that FBCs, indeed, hmit the de- 
coherence and relaxation, one can explain [l6| experi- 
mentally observed [l3| 1// noise only provided that the 
number of FBCs is large. Thus it looks that the above 
considerations underestimate the efficiency of this source 
of decoherence. 

In this letter we show that allowing for the Coulomb 
interaction between the trapped charge and bath elec- 
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trons results in a dramatic increase of the effective hy- 
bridization rate 7 and thus in the upward revision of the 
estimate ([T]). We wih present fuU analytic results that 
show how the interaction enhances decoherence. 

We consider the model [l^ [ill of the charge qubit 
indirectly coupled to the thermal bath via the coupling 
^ ^ X^i a ^i^icr^ia FBCs hybridized with the bath: 
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Here wq = Ec{N — Ng) is the energy gap between the 
qubit levels, Ti are the Pauli matrices in the space of 
qubit states, E is the control energy (that includes the 
Josephson energy in case of the JCQ); , d^^ are the 
creation and annihilation operators for an electron with 
spin a =tj i the j"^ trap and Vi is the electrostatic 
coupling of the qubit to the trapped electron. 

We assume that the bath, described by i/e, is a nor- 
mal metal - the situation which is more relevant for the 
DQD charge qubit (for the JCQ this mechanism might 
be relevant as well since in the absence of screening spa- 
tially remote trapped electrons hybridized with metallic 
gates may contribute to decoherence). In this case 



Hn — Hn 
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Here H„ 



J2p a ^p^p.aCp,a the Hamiltouian of the 



metallic bath, c^, ^ , Cp ^ creation and annihilation oper- 
ators for bath electrons. The energy Sp = /2m — ep 
is counted from the chemical potential (and so are be- 
low). The impurity Hamiltouian describes random levels 
which can trap electrons. 



E 



H- 
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where hirr — d\^d^^ and n. 



ni-\-\-nii. The Hubbard on- 
site repulsion C//^ is assumed to be very large so that two 
electrons cannot reside on the same site. The relevant 
charge fluctuations are due to switchings between empty 
and single-occupied impurity states. In what follows we 
disregard double-occupied states and omit irrelevant spin 
indices, absorbing the spin degeneracy of the bulk elec- 
trons into the density of states = 1/V5. 

The d-electrons are coupled to the bath via direct hy- 
bridization as well as via the local Coulomb interaction: 
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Here ti is the tunneling amplitude from the resonant level 
at the i^^ trap into the bulk, pi is the bath electron den- 
sity operator at the i^^ site and the interaction Ui is pos- 
itive for an acceptor and negative for a donor level. 
The model of E qs. ( Et- E l) was previously considered 
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tions Ui in Eq, 
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M\2M,\21,\M without the interac- 
We will extend the non-perturbative 



procedure developed in [l8| to evaluate the contribution 
to decoherence of the charge qubit due to FBCs interact- 
ing with the bath electrons. 

Let the full density matrix of the system initially be 
separable, i.e. /3b+q(0) = p{0) (8) Pb, where p{t) is the 
reduced density matrix of the qubit and pB = Z^^e'^^^ 
is the equilibrium density matrix of the bath. Then one 
can formally represent p{t) as 



m 



n(t) 



Pi2{0)e-'^'>*V{t) 



P2i{0)e'^«'V*{t){t) l-n{t) 



(6) 



In the long-time limit, both diagonal and off-diagonal 
elements of p(t) get exponentially close to their equilib- 
rium values. The corresponding decay rates are usually 
referred to as the relaxation, Fi, and decoherence, F2, 
rates. Below we discuss "pure dephasing" , i.e. let £' = 
in Eq. ([2|). This is justified for the charge qubit during 
the most of the operational cycle 15]. We find the time 
evolution of 'D{t) as 



AHB + V)t -i(HB-V)t 
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where {A) = Tr(psi). As F2 = - limf^^o t-^ln\D{t)\, 
we need to calculate the asymptotical value of 2?(t) at 
i — > 00. To this end, we represent 2?(t) in Eq. ([7]) as a 
functional integral over the Keldysh double-time contour, 
which is natural as Eq. ([7|) contains time-ordered and 
anti-time-ordered exponentials. 

Note now that Eqs. ([3|)-([5|) determine a well known 
problem of the interacting resonant level (IRL) which can 
be mapped [2^ onto the single channel Kondo problem. 
This suggests that one can account for the interaction 
(the second term in Eq. (I5|)) by using the renormaliza- 
tion group approach. We will show that this leads only 
to substituting the renormalized (energy-dependent) hy- 
bridization rate for the bare one into Eq. Ijllap below. 

In the absence of the interaction term in Eq. ([5]) one 
performs an exact functional integration over the fields 
describing electrons in the bath. This results in the ef- 
fective action in terms of the fields di{t), di{t) which 
describe trapped electrons 

S[d\d] = Yl Jdtdt'dl{t)G^\t,t')d,{t') . (8) 

The time integrations are performed along the standard 
Keldysh contours. The Keldysh Green function Gij{t,t') 
of the FBC is defined by the Dyson equation resulted 
from the integration over the bath electron fields: 

G^\t,t') = S{t,t')d,j[idt-e,+v,{t,t')] - %j{t,t'), 
t,j{t,t')^Ut*g{n-r'j-t-t'). (9) 

Here Vi{t, t') = ±d{t)9{t — t')vi with sign on the upper 
and ' — ' sign on the lower branch of the Keldysh contour. 
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and g is the bare Keldysh Green function of the con- 
duction electrons: the Fourier transforms of its retarded 
component is g^{e,p) = {e — Ep + iO)^^. Typically, FBCs 
are separated by a distance far exceeding Ap- Then the 
off-diagonal contributions to is suppressed upon a 
spatial integration. Thus = GiSij with Gf (e) = 
{e - ei+ij,/2)-^ and Gf(e) = 2ImGf ((r)tanh(£/2r). 
The width 7^ = ReT^u = 7r|iip/i5 of the i^^ trap level 
plays the role of the recharging switching rate. 

Since contributions of different FBCs are uncorrelated, 
we consider a single fluctuator with the switching rate 70. 
Then the functional averaging in Eq. ([7]) with the action 
(111) results in 



To 



-Re lim Tr ln(l + wG). 

t — >oo 



After expanding this in powers of vG , one finds that 
the factor cancels while the remaining time-integrals 
in each order have the convolution structure (3]. This 
allows to re-exponentiate the expansion after perform- 
ing the Fourier transform. Taking then the trace in the 
Keldysh space, we arrive at 

-t-oc 

T2^-^j ^ln[l + wG^(e) -i;G^(e)«G^(e)] . (10) 



Substituting G^-^ and G^ into Eq. ^ we find T2 as a 
function of T and 7(e): 



r2 = -/ ^ln{l-4/(e)[l-/(e)]sinV(e)}, (Ha) 



where /(e) = (1 + cF/'^) ^ is the Fermi factor and 



^p{e) = arctan 



2v"f{e) 



(e - eo)^ + 7^(e) - v'^ 
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Note that the interaction in Eq. ([5]) entered r2 only 
through substituting the renormalized energy-dependent 
switching rate 7(e) for the bare value 70. The RG equa- 
tions for the IRL model, Eqs. (I3|)~(l5|), have been ob- 
tained by mapping the tunneling amplitude to onto the 
transverse coupling in the Kondo model (see, e.g., 
[i^ilHill^l). These RG equations are equivalent to those 
originally derived for the Kondo model within the ID 
Coulomb gas picture (soj . It is convenient to write them 
in terms of the effective charges Y and (5, 
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where 5 = 2 arctan (7ri'o?7/2), Y = (7/e)cos^^ and e is 
the running energy cutoff. The upper cutoff is the band- 
width D of the bath electrons. As typically 70 <C the 
bare value Yq ~ 7o/-D ^ 1 so that one can neglect the 
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FIG. 1: Decoherence rate as a function of coupling strength, 
g = v/"f, for T = O.I7, D — IO7. The dashed line shows the 
non-interacting case {a — 0), while the solid line describes de- 
coherence with allowance for the interaction (a — 0.5) which 
turns out to be substantial in a much wider parametric range. 



renormalization of 5. Substituting S = Sq, we rewrite the 
RG equations as 



din 7 
dine 



2^ 

TT 
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Equation p2)) yields 7(e) = ^o{D/e)". The energy de- 
pendence saturates for small e at e — max(7, T,ea). The 
parameter a can be positive or negative, depending on 
the interaction sign. For Uq > the linear in S contri- 
bution to a (due to the Mahan exciton) is dominant 31 1 
and exceeds the quadratic contribution (due to the An- 
derson orthogonality catastrophe) so that a > 0. Below 
we take into account only the levels with positive a. 

According to Eq. (fTT|) all traps with 7(e) > Eq con- 
tribute effectively to Since 7(e) » 70 for a > 0, 
the fraction of effective FBCs dramatically increases in 
comparison to the non-interacting case. 

At a given T each FBC is either in the high- 
temperature, T > 7(e), or in the low-temperature, 
T < 7(e), regime, depending on the relation between 
T and To = 7o(-D/7o)^. As eo ^ 7o the effective 
number, Eq. ([1]), of "high-T" FBCs can be estimated 
as (T'/D)^"'""7Vgeom, making their presence rather im- 
probable - but not impossible. The switching rate for 
such FBCs is T-dependent and saturated at 7(T) = 
-foiD/T)". Most (if not all) FBCs are in the low-T 
regime, T < 7(e), i.e. characterized by To exceeding 
T. In this case the switching rate 7(e) saturates at 

/ 7-, / \a/{l+a) 

7 = 70 {D/j) 

Since at low T the product of Fermi factors in Eq. (|lla[) 
is a peak around zero (i.e. the Fermi energy) of width T, 
we approximate 7(e) by the £-independent expression, 



7 « 7(e = 7o) 



_ / 7o(-D/7o)°+^ 



7 £0 



lo{D/eo) , 7^eo 



(13) 
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Thus the contribution of the low-T FBC into the deco- 
herence rate is 



r^°"(r) = -arctan^ 

TT 



2W7 



4+T 



(14) 



Figure 1 illustrates this behavior and shows that the in- 
teraction dramatically increases the energy where FBCs 
effectively contribute to For Eq ^ 7 the decoherence 
rate decays with Eq faster than in the non-interacting 
case: r2 oc e^'*"^". This means that, as was noticed 
above, only FBCs with eo ^ 7 cause decoherence. Using 
7 from the upper line in Eq. ()13|) we find the effective 
number of FBCs: 



off ' 



(15) 



This is considerably higher than that in the absence of 
the interaction, Eq. ([1]). 

We have already noticed that it is statistically improb- 
able to find the FBCs in the high-T regime. Had they 
existed, their contribution would be also enhanced by 
the interaction as compared to the expressions obtained 
classically 1^ 3, namely 



(16) 



where 7(T) = 70 (I?/T)". Eq. (fTB|) predicts a non- 
analytic in T decoherence rate with a well pronounced 
peak at T = Dg'^/"'. 

In conclusion, we have demonstrated that the fluctu- 
ating background charges - namely, electrons trapped 
in the insulating layer close to the surface of metal- 
lic gates/electrodes and thus hybridized and interacting 
with the conduction electrons - cause significant deco- 
herence of the charge qubit. The interaction leads to 
a drastic increase of the efflciency of this mechanism of 
decoherence and to a parametric enhancement of the re- 
sulting contribution, as well as of the number of the effec- 
tive FBCs given by Eq. (jisp . This enhancement proves 
the assumption that the FBCs act as the most dangerous 
source of decoherence for the charge qubit. 
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